Soft elasticity of RNA gels and negative Poisson ratio 
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We propose a model for the elastic properties of RNA gels. The model predicts anomalous elastic 
properties in the form of a negative Poisson ratio and shape instabilities. The anomalous elasticity 
is generated by the non-Gaussian force-deformation relation of single-stranded RNA. The effect is 
greatly magnified by broken rotational symmetry produced by double-stranded sequences and the 
concomitant soft modes of uniaxial elastomers. 
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I. INTRODUCTION 

The study of the viscoelastic properties of networks 
of flexible, synthetic polymers has for many years been 
a central topic of polymer science. The classical Flory 
theory for the elasticity of rubber and of gels treats 
these systems as networks of nodes linked by highly flex- 
ible chains 0, 0. Scaling relations for the viscoelas- 
tic moduli that result from this model have been well 
connrmed[j|- The study of networks of biopolymers has 
provided a fresh impetus to the held. Gels of semi- 
flexible biopolymers, like F-Actin, were shown to obey 
novel scaling relations The focus of the present pa- 
per is on the elasticity of a different biopolymer system, 
namely a network or gel of RNA chains. The folding 
of smaller RNA molecules has been already extensively 
discussed in the molecular biology literature in the con- 
text of Ribozymes[5|, but extended RNA gels have not 
received much attention. The genome of single-stranded 
(ss) RNA viruses for instance may form a very promising 
small scale realization of an RNA gel, as discussed in the 
conclusion. R is the purpose of this paper to present a 
simple model for RNA gels that indicates that such gels 
should have rather unusual elastic properties that will 
distinguish them not only from Flory-type gels but also 
from gels of semi-flexible biopolymers like F-Actin. 

A single-stranded ("ss") RNA chain can be folded 
first into a "secondary" structure that consists of the 
pattern of optimal pairing of the bases of the chain 0. 
This secondary structure is represented as a planar, 
branched graph of duplexed double-stranded ("ds") se- 
quences linked by "bubbles" and "stem-loops" composed 
of unpaired bases. A three-dimensional "tertiary" struc- 
ture is obtained if one also allows complementary pairing 
between the bases of different bubbles and stem-loops 
of the secondary structure. The model we will study as- 
sumes a highly simplified tertiary network topology com- 
posed of an array of rigid rods (the ds complementary 
sequences) that are linked by flexible chains (the ss se- 
quences). Two flexible chains emerge from either of the 
two ends of each rod (see Fig. Models of this type 
have in fact been used to describe the folding kinetics 
of Ribozymes |6| . A key ingredient of our model is that 
these flexible chains are not assumed to have the elastic 
properties of either a Gaussian chain or a semi-flexible 



worm-like chain. Instead, we will examine the elasticity 
of the network for a general class of interaction poten- 
tials. Specifically, we will consider chains that obey the 
force-extension curve of ss DNA as measured by single- 
molecule micromechanics. 

The particular network of rods and springs that we 
propose to investigate as a representation of complcxcd 
RNA is displayed in Fig. 2] The rods are centered on 




k — 2, 



FIG. 1: The network of rods and springs that forms the basis 
of the negative Poisson ratio system. As indicated in the 
figure, the length of the sides of the underlying altered cubic 
lattice is a and the length of the rods is 21. 



an altered square lattice. The angle, 4> that an edge of 
this lattice makes with respect to the vertical, as shown 
in Fig. parameterizes the overall characteristics of the 
lattice. If cf>, which ranges between and 7r/2, is equal to 
7r/4 the lattice is square. At either of the limits of (j>, the 
lattice has collapsed onto itself. The edges of the unit 
cells of the lattice, which do not represent any physical 
quantity, are represented by grey lines in the figure. The 
rods are thick solid lines, and the springs are depicted 
as dashed lines connecting the ends of neighboring rods. 
As shown in the figure, the length of the edges of the 
underlying lattice is a, and the length of the rods is 21. 
The model belongs to a class of systems, namely uni- 
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axial and biaxial elastomers, that are known to have un- 
usual elastic properties. The internal rotational degrees 
of freedom produce what is known as "soft" elasticity in 
the form of large shape changes under applied fields as 
well as a vanishing of the Poisson ratio The anoma- 
lous elasticity of networks with broken rotational symme- 
try is an unavoidable and fundamental feature according 
to a theorem by Golubovic and Lubensky 8] . One of the 
aims of the present paper is to examine the consequences 
of this theorem for a concrete model. Specifically, we will 
explore precisely what physical properties of the interac- 
tion potential characterizes anomalous soft elasticity. We 
will argue in particular that control of the soft elasticity 
can be achieved by altering the physical properties of the 
non-Gaussian springs and that this can be achieved in 
the context of RNA gels. 

The proposed model can be viewed as a modified ver- 
sion of a hexagonal lattice of hinged rods. It is actually 
known in the material science literature that systems that 
can be represnted by models of this type exhibit unusual 
elasticity. A conventional hexagonal lattice, shown at the 
top Fig. [21 will expand vertically when compressed hor- 
izontally if the rods are rigid and freely hinged. On the 




FIG. 2: Fwo types of hexagonal networks of hinged rigid rods. 
The network shown at the top of the figure exhibits a positive 
Poisson ratio, in that it expands vertically when compressed 
horizontally. The network at the bottom contracts vertically 
when compressed horizontally and thus has a negative Poisson 
ratio. 



other hand, the network at the bottom of the figure — 
consisting of non-convex hexagonal units — contracts ver- 
tically when compressed horizontally. The network is 
thus characterized by a negative Poisson ratio. The con- 
ditions for constructing a network of hinged rods that 
exhibits a negative Poisson ratio can be elegantly demon- 
strated by the popular Hoberman Sphere [!j , which main- 
tains its spherical shape as it expands and collapses [10). 
The generic term for a material with a negative Pois- 
son ratio is "auxetic." A negative Poisson ratio is known 
to be exhibited by Iron Pyrites [llj, self-avoiding, fixed- 



connectivity membranes [12|, monocrystalline Zinc |13| . 
Carbon Nitride [l4| , polyethylene foams 0] , two dimen- 
sional mesh-like s yste ms fl6l| , structures composed of ro- 
tating rigid units |l7| and in network-embedded compos- 
ites [l8j. Biologically, a negative Poisson ratio has been 
found to be exhibited by both skin and bone poT j. 

II. POISSON RATIO OF A TWO 
DIMENSIONAL NETWORK OF RODS AND 
SPRINGS 

We will start by exploring the properties of the model 
shown in Fig. as a two-dimensional network. We will 
represent a distortion of the lattice in terms of a strain 
tensor e , so that the displacement Af^ = (Axi, Ay^), of 
a lattice vertex originally at the location fj = (xi,yi) is 
given by 

AXj — txx-E'i ~h txyVi (^--0 

and similarly for Az^. We will also allow the rods to 
rotate in the plane by the angle 0. 

Consider, now, the triad of rods with connecting 
springs shown in Fig. [3] Note that the underlying lattice 




FIG. 3: A triad of neighboring rods, connected in this figure 
by two springs connected to the top end of the lower cen- 
tral rod. The rods and the lattice are shown in an as-yet 
undistorted state. The lattice is compressed horizontally as 
compared to the one shown in Fig. Q 

is more compressed along the horizontal direction than 
the lattice shown in Fig. ^ As will be demonstrated 
below, the degree of alteration of the square lattice is 
controlled by the requirement that the rod and spring 
network can be stabilized by an osmotic pressure. The 
distance between the upper end of the lower rod and the 
lower end of the rod on the right is equal to 

d\ = (ae xy cos (j) — 21 sin# + a sin (f> + ae xx sm(j>) 2 
+ (—21 cos 9 + a cos <j) + ae yy cos <\> + ae yx sin 4>) 2 

(2.2) 



3 



and the distance between the upper end of the lower rod 
and the lower end of the rod to the left is given by 



d 2 = (a£a;j,COS( 



21 sin 6 



a sin m — ae xx sin < 



+ (— 2/cos6' + acoscj) + ae yy cos< 



(2.3) 



The angle 9 in 112.211 and l|2.3|l is the angle to which all the 
rods rotate under the influence of the uniform distortion 
of the lattice as parameterized by the constant strain 
tensor e . 



A. Non-Gaussian behavior of the springs 



all higher order coefficients in the expansion of the inter- 
action in terms of the variable x. On the other hand, for 
a freely jointed chain or a worm-like chain V2 > 

Under the assumption that the interaction mediated 
by the springs is always attractive, the angle, 9 at which 
the rods tilt must adjust in such a way as to minimize 
the sum of the two distances, d\ and d%. Taking the 
derivative of d\ + g?2 with respect to 9 and constructing 
an extremum equation, we end up with the relationship 



8al cos </> [(e 



1) sin# — e xy cos( 



= 



(2.9) 



The strain tensor is assumed to be a small quantity. This 
means that the solution to the equation above is, to an 
accuracy sufficient for our purposes 



In order to facilitate numerical calculations, we will 
parameterize the interactions mediated by the springs in 
terms of the square of the distance between the spring 
end-points. That is, we express the energy of the inter- 
action in terms of the variable x, where 



dl 



(2.4) 



The quantity d is the distance between end points, and 
do is the distance between end-points in the equilibrium 
state. Next, let V{d) be the interaction potential that 
represents the free energy of an ss sequence connecting 
two points a distance d apart. We will expand V(d) in a 
Taylor series around the equilibrium state x = 0: 



V ( xldl + x 



= v{do) + ^~ v ' {do) + If [ ~ v ' {do) + d ° v "^ 



H 

vq + V\x + —x 



(2.5) 



It is important to note that even if the interaction, V(d), 
is convex upward, so that V"(d) > 0, the second order 
coefficient, V2, in the expansion above can be negative. 
Consider for instance the case of a simple power-law in- 
teraction energy of the form 



V(d) = CJP 



(2.6) 



with C a positive coefficient and p a power greater than 
one. Then, 



and 



V'(do) = Cpd^ 1 (2.7) 



d V"(d ) - V'{d ) = Cp( P - 2)d p ~ 1 (2.8) 



If p < 2, then, assuming that the coefficient C is greater 
than zero, the second order coefficient in l|2.5(l . will 
be negative. In the case of Gaussian chains with Hooke- 
law-type harmonic interaction, V2 is equal to zero, as are 



1 rtj 



(2.10) 



The total energy associated with the interactions me- 
diated by the springs shown in Fig. [21 is then given by 



E 



v Q + vx ((dl - dl) + (4 - d 2 Q )) 

+ v i((di-dif+(di-di) 2 



(2.11) 



with 9 in and (JOJ) as given by l|2.10|l . 



B. Determination of the angle cj>, expansion of the 
energy in strain coordinates and calculation of the 
Poisson ratio 

The next step is to expand the resulting expression to 
second order in the strain tensor. At zeroth order, we are 
left with the coefficient vq. To first order in e we have 
the following contribution to the energy 



Aav\ (ae xx sin 



i^)-2l)) (2.12) 



The appearance of terms linear in the strain energy in- 
dicates that the system has a tendency to spontaneously 
deform. First, consider the case of a uniform contraction. 
In terms of the strain tensor, this yields an energy going 
as e xx + E yy . This tendency will be countered by an os- 
motic pressure II that increases as the monomer concen- 
tration grows under contraction. In order that an osmotic 
pressure can completely balance the first order energy for 
d = do (the equilibrium state) in l|2.12|l . it is necessary 
that the coefficients of e xx and e yy in that expression be 
equal to each other. This requirement translates into an 
equation for the lattice angle <f> that is satisfied when 



arccos 




(2.13) 



If the coefficients of e xx and e yy are unequal, then the sys- 
tem undergoes a spontaneous shear tranformation. Two 
limits of the above expression are noteworthy. When 
I = 0, so the rods are infinitesimal in extent, then 4> as 
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given by i|2.13[) is equal to tt/4, consistent with an under- 
lying square lattice. On the other hand, when I = a/2, 
we find a 4> from (|2.13(1 that is equal to zero, which im- 
plies a horizontal collapse of the complex. In this limit, 
the rods are exactly long enough that their tips touch in 
the event of such a collapse of the lattice. Henceforth we 
will assume that I lies in the range between and a/2. 

Assuming cancellation of the linear term, the energy 
is now, at lowest non-trivial order, quadratic in the ele- 
ments of e . This contribution to the quadratic energy is 
of the form 



(2.14) 



In light of translational invariance, this energy will be the 
same for the energy supplied by the two springs attached 
to the upper tip of every rod in the complex, under the 
assumption of a uniform strain. We have thus effectively 
calculated the energy of interaction in the entire complex. 

The strain induced by an externally applied stress, cr, 
is the solution to the set of equations 



d 



(2.15) 



We will take the strain to be exerted in the y direction, 
so <jyy is the only non-zero element of the strain tensor. 
After solving 12.15fl for the elements of the strain tensor 
we can extract the Poisson ratio, /i, via the relationship 



fj, 



(2.16) 



!)!) 



Following a straightforward calculation we obtain the ex- 
plicit result 

a 2 (a 2 -4Z> 2 



(a 2 - I 2 + aly/2 + (l/a) 2 )^ + a 2 (a 2 



4Z> 2 
(2.17) 



Poisson ratio when there are no rods 



lattice, provided — vi/3a 2 < v 2 < i^i/a 2 , where we as- 
sume that vi > 0. For v 2 in this range, the Poisson ratio 
varies from -1/2 to 1/2. It follows that a moderately 
negative Poisson ratio is, in fact, possible in a lattice 
consisting entirely of non-Gaussian springs. Note that 
this is not related to the existence of internal rotational 
degrees of freedom. 



D. Poisson ratio at general I 

We now set I — 0.45a. This leads to an equilibrium 
state that is compressed horizontally as in Fig. [3] In 
Fig. 0]we show /i as a function of v 2 . 
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FIG. 4: The two Poisson ratio of the two-dimensional lattice 
according to (12. 1711 . 



We have again checked the stability by computing the 
determinant of the four-by four matrix with elements 
ctij^ki- Figure shows a typical plot of the determi- 
nant as a function of v 2 . The determinant vanishes at 



Det(SE) 




v, 



It is useful to consider certain limiting cases of Eq. 
(|2.17l) . First, consider the case I = 0, which corresponds 
to a network in which the rods are replaced by point 
contacts. If we set I = in the expression l|2.17[) . the 
expression for the Poisson ratio further simplifies to 



a 2 v 2 /vi 
1 + a 2 u 2 /v\ 



(2.18) 



which predicts a negative Poisson ratio, if the ratio v 2 / v\ 
is negative. To check for stability of the lattice one can 
determine the eigenvalues of the energy matrix a . When 
I = 0, this four-by-four matrix is degenerate, and there 
are two distinct eigenvalues: a 2 Vi—a i v 2 and a 2 i^i+3a 4 ^2- 
Both eigenvalues are positive, corresponding to a stable 



FIG. 5: The determinant of the elastic energy matrix a when 
I = 0.45a, and v\ and a have been set equal to one. 



v 2 = —6.232, corresponding to the emergence of a nega- 
tive eigenvalue of the energy matrix a and the develop- 
ment of a mechanical instability in the complex |2l| . At 
that threshold value of v 2l we find fi = —4.21226, so a 
substantial negative Poisson ratio is indeed possible close 
to a mechanical instability. Figure shows the determi- 
nant for an extended range of the interaction coefficient 
v 2 . As illustrated in that figure, the range of stability 
of the energy matrix is bounded from above as well as 
below as a function of v 2 . 



5 



Det(a) 




FIG. 6: The determinant of the energy matrix a as a function 
of V2, with I — 0.45a and a and v\ set equal to one. The range 
of stability is indicated by the heavy line on the horizontal 
axis. 



If we compare these results with those in the case 
! = we conclude that the broken rotational symme- 
try for I has greatly amplified the negative Poisson 
ratio produced by a negative value of vi. There is in- 
deed no intrinsic physical bound on the Poisson ratio in 
an anisotropic solid [22] • The amplification effect is most 
dramatic when vi approaches the threshold of the me- 
chanical instability, at V2 = —6.232. Note though that 
according to Eq. (|2.17|1 the Poisson ratio formally is zero 
even at this critical point for a Gaussian network and any 
other system with ^2 = 0. It should be recalled here that 
the Poisson ratio is zero in conventional nematic elas- 
tomers Q. One can compare the effect of a negative v-i 
to that of a small magnetic field applied to a magnet as 
one approaches the Curie temperature. The divergence 
of the susceptibility amplifies the effect of the applied 
field. In the next section we will extend the notion of 
the sort of complex we have been discussing to a three 
dimensional system. 



III. THREE DIMENSIONAL NETWORK OF 
RODS AND SPRINGS EMBODYING A 
NEGATIVE POISSON RATIO 

A. Preliminaries: description of parameters 

The way in which the three-dimensional network 
is constructed will be a generalization of the two- 
dimensional case. Imagine a sheared version of a cubic 
lattice rotated so that the z-axis is along a diagonal. A 
rod is placed at each vertex of this lattice, and each end 
of this rod is connected to the ends of the three rods that 
are closest to it by a spring having a non-linear force- 
extension relation. Figure depicts a portion of this lat- 
tice of rods, shown as short solid lines. The springs are 
represented by dashed lines and the edges of the lattice 
shown for illustrative purposes only — are represented by 
grey lines. The z axis is along the principal diagonal, 
shown as a grey dashed and dotted line in the figure, and 
the x axis is aligned with the projection in the x-y plane 
of an edge of the cube. Three neighboring vertices of 




FIG. 7: Portion of the cubic lattice, shown unsheared here, 
with the rods at the vertices. The "springs" connecting the 
vertices are shown as dashed lines. The edges of the cubic 
lattice in which the rods are embedded are shown as grey 
lines. Note that these edges do not represent actual physical 
structures. Only those springs that connect two rods shown 
in the figure are depicted. Finally, a principle diagonal, lying 
along the z axis, is shown as a dashed and dotted line. 



the cube are indexed as shown in Fig. |SJ The organiza- 



21 




FIG. 8: Three vertices of the cubic lattice with associated rods 
and springs attached to them. The numbers are the indices 
of the vertices. As in Fig [7] the springs are shown as dashed 
lines. 



tion illustrated there ought to be understood as a kind 
of "averaged" version of a structure in which each rod 
terminates in two springs, as shown in Fig. ^ 

The rod with index is assumed to be located at 
the origin. The coordinates of the corners on which the 



6 



neighboring rods sit are given by 



Xi 


= a cos ip 


(3.1) 


Vi 


= 


(3.2) 


Zl 


= a sin ip 


(3.3) 


X2 


= A cos ip cos(2tt / 3) 


(3.4) 


1)2 


= a cos ip sin(27r/3) 


(3.5) 


Zl 


= a sin ip 


(3.6) 


X3 


— a cos ^008(4^/3) 


(3.7) 


2/3 


= a cos ^ sin (47r/3) 


(3.8) 


Z3 


= a sin ip 


(3.9) 



The quantity a in Eqs. H3.1I) - (|3.9|) is the length of a cube 
edge. The quantity ip is the angle that the edges of the 
cube in Fig. |H1 make with respect to the x-y plane. In 
the case of an unsheared cubic lattice, this angle is equal 
to arctan(l/v / 2). Note that when ip = tt/2, the edges 
connecting the site at the origin to the three neighboring 
ones are all vertical, so that the neighboring sites all lie 
at exactly the same location. Each edge in the cubic lat- 
tice has length a and the length of the rods is 21, again as 
shown in Fig. [S] Initially, the rods point along the z-axis. 
However, we will assume that when the lattice shears the 
rods change orientation. The direction in which the rods 
point will be described in terms of the standard spherical 
angles 9 and (p. See Fig. E|for the representation of the 
angles <p and 8. The next step is to allow the lattice to 




FIG. 9: The angles cp and 6. All three edges shown in the 
figure have an orientation with respect to the horizontal plane. 



distort. As in the previous section we represent the dis- 
tortion in terms of a strain tensor e , so that placing one 
of the vertices of the lattice at the origin, the displace- 
ment Afi = (Axi, Ay.i, Azi), of any other lattice vertex 
originally a the location fi — (xi, yi, zi) is given by 



AXi — txxXi ~\~ txyVi ~\~ ^xzZi 

and similarly for Ayi and Azi. 



(3.10) 



B. Determination of the angles ip, (p and 9. 

First, we calculate the quantities (xi, yi, Zi) for each of 
the springs shown in Fig. |SJ Then, we expand the energy 
to first and second order in the (xi, yi, Zi). The resulting 



expression is then expanded once again, now to second 
order in both the strain tensor e and the rotation angle 8. 
As both quantities are expected to be small, we restrict 
our consideration to terms that are at most quadratic in 
either one or both of them. The next step is to determine 
the angle ip, (shown in Fig. [5J by which the lattice shears. 
This angle is again controlled by the requirement that 
the lattice, in the absence of any imposed stress, can be 
stabilized by an osmotic pressure II, that couples only to 
the total volume of the lattice, or in other words to the 
combination Tr e = e xx + e yy + e zz . 

If we expand the energy of the lattice to first order in e 
and zeroth order in 8, we obtain the following expression 

Sa 2 {£xx + £yy) cos 2 ip + 6ae zz sintp(asmip — 21) (3.11) 

The requirement that this reduce to a function of the 
trace of the strain tensor only leads to a quadratic equa- 
tion for smip. The solution to this equation is 



iP 



21 + V30 2 



4/ 2 



3« 



*o(//a) 



(3.12) 



It is useful to consider the two limits of the above expres- 
sion. When I — ► 0, the argument of the arcsin reduces 
to 1/ V3, which is consistent with an arctangent equal to 
1/V2 and a cubic lattice. On the other hand, as I — > a/2, 
the angle ip goes to 7r/2, which means that the lattice col- 
lapses onto itself. In the latter limit, the rods are long 
enough that they touch end-to-end when that collapse 
takes place. 

Inserting the value of ip given by the right hand side 
of l|3~T2")l into the expression H3.11|l . we find for the de- 
pendence of the energy on uniform dilations of the lattice 



2 , 

3 y^xz 



3a 2 -4Z 2 -2ly / 3a 1 



so the Equation of State of the system is 
2 



n« 



3a 2 -4^ 2 -2ZV30 2 



Af 



(3.13) 



(3.14) 



where 



= a 3 (3\/3/2)cos 2 ?/>smV> 

(3a 2 - 8Z 2 )V3a 2 + 4Z 2 - 16Z 3 



(3.15) 



is the volume per rod in the lattice. As / — * a/2, the 
expression multiplying the trace of the strain tensor van- 
ishes, and the linear dependence of the energy on uniform 
dilations of the lattice is lost. The osmotic pressure, n 
approaches the limiting value 2v\l\pha. 

In light of this initial adjustment of the lattice, the 
next step is to determine the extent to which the lattice 
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responds to an external stress. The response manifests 
itself in three quantities: the strain tensor e , the rod 
tilt angle 9 and the azimuthal angle, tp, of the rods. As 
it turns out, the azimuthal angle appears only in the 
term that is first order in 9, and the form of that term 
is A cos 4> + B sin cj), where A and B are linear functions 
of the strain tensor and general functions of the angle tp. 
This expression is minimized when (f> — arctan^l/i? and 
is equal in this case to — V A 2 + B 2 . The second-order- 
m-9 contribution to the energy of the deformed lattice is 
equal to 



Gal 



f = D 2 t 



2alv 2 (cos* (V°)) + v i sin* (i/ a ) 

(3.16) 

If we denote by Do the portion of the deformation enerrgy 
that is zeroth order in 9 and second order in the strain 
tensor, the net free energy as a function of the polar angle 
is 



D - 9^ A 2 + B 2 



D,9 2 



(3.17) 



Minimizing l|3.17|l with respect to 9, we end up with our 
final result for the dependence of the energy on the strain 
tensor e in the uniformly strained lattice. In terms of 
the quantities defined above, the energy has the form 



D - 



A 2 +B 2 



4L> 2 



^ OiijMtijtkl 
i,j,k,l 



(3.18) 



As indicated in (|3.18(l , the energy of the uniformly dis- 
torted lattice with the angles ip, <f> an d 9 relaxed to the 
values dictated by energy minimization and lattice sta- 
bilization is purely quadratic in the strain tensor. The 
tensorial quantity cx-ij^i is symmetric with respect to in- 
terchange of the index pairs, ij and kl. The response of 
the system of rods and springs to an externally generated 
stress tensor, cr, follows from the solution to the set of 
linear equations 



(3.19) 



Given the overall anisotropy of the lattice, symmetry ar- 
guments yield only restricted information. If we orient 
the z axis along the direction defined by the rods in the 
unstressed lattice, then symmetry arguments tell us that 
a stress entirely in the z direction results in equal val- 
ues of the strain tensor components e 



and e yy , while 



C. Poisson ratios for / = 

Once again we consider the case of a lattice consisting 
of springs only, in which the rods are replaced by point 



connections. For I — 0, the Poisson ratio of interest be- 
comes 



4(a 2 is 2 /vi) 



and the determinant of the energy matrix is 



<l.<t n=:»L2a 1 Vfl + — " 



(3.20) 



(3.21) 



This means the energy matrix is associated with a stable 
equilibrium state as long as a 2 v 2 jv\ > —1/2. FigurelTUlis 
a plot of the Poisson ratio as given by l|3.20|l in the range 
of stability of the energy. The limits of this quantity are 




a 2 v 2 A^ 



FIG. 10: The Poisson ratio as a function of a 2 V2/vi when 
I — 0. The asymptote of one is also indicated. 



1 as a v<xjv\ — > oo and -2 as a 2 v 2 jv\ — > —1/2, which 
are the limits of the Poisson ratio in an isotropic solid. 
This means that, as in two dimensions, broken rotational 
symmetry is not essential for the negative Poisson ratio. 



D. Poisson ratios for / 7^ 

To explore the properties of the Poisson ratio for I ^ 0, 
we will focus on the exponent p in the power law interac- 
tion l|2.6(l . Consider the graph shown in Fig. ^2 which 
is of the Poisson ratio, defined as 



(Jl : 



c yy 



(3.22) 



in the case of a stress entirely in the z direction (er^- = 
except for a zz ). The length, Z, of the rod has been 
set equal to 0.4a. Note that the Poisson ratio passes 
through zero as the power law passes through the power 
associated with the Gaussian spring. When p < 2, the 
Poisson ratio is negative, which means that the lattice 
system resists shearing. 

FigurelTUlis a graph of the Poisson ratio when I = 0.45a 
as a function of p. The range of stability, i.e. of a pos- 
itive determinant for a, is for p > 1, with p = 1 the 
location of the onset of instability. Recall that in the two 
dimensional case, fj, approaches a finite negative value at 



8 



1 

s 

0.2 


■ 

1 


1.6 1.6 j*gf& 
_ .2 


2-2 2.4 

- 


/ -0.4 




/ -0.6 





p 



FIG. 11: Graph of the Poisson ratio versus the exponent, p, 
of the power law in 12. 611 . 




FIG. 12: The Poisson ratio, plotted against the exponent p 
m [THfll . Here p ranges from zero to three. As noted in the 
text, a value of p less than one is inconsistent with mechanical 
stability. 



the point of instability. From Fig. ^| we see that the 
Poisson ratio diverges to minus infinity at the instabil- 
ity threshold. The amplification effect of the rotational 
degree of freedom in producing a large negative Poisson 
ratio is thus much more pronounced in three dimensions. 
By "tuning" the physical properties of the non-Gaussian 
springs it is possible to dramatically alter the elastic re- 
sponse of the network. 



IV. RNA IN NETWORKS AND A NEGATIVE 
POISSON RATIO 

In this conclusion, we will apply the results of the last 
section to RNA networks and discuss the consequences 
in terms of soft elasticity and negative Poisson ratios.. 
In Sections [H] and 1 1 1 1 1 we learned that functions as 
a control parameter for elasticity. If is positive, one 
obtains conventional elasticity. As noted in the Intro- 
duction, gels of flexible and semi-flexible polymers are 
conventionally modeled as networks of Gaussian chains, 
worm-like chains or freely jointed chains, all of which have 



V2 > 0. If V2 is negative, we predict a range of anoma- 
lous elasticity, with negative Poisson ratios terminating 
in a mechanical instability at a critical value for at 
which point ji diverges to minus infinity. Is it realistic 
for a biopolymer network to have a negative The 
Vi parameter can in principle be determined from the 
force-extension curves of the nonlinear springs. For an 
RNA network of the type shown in Fig. ^ this would be 
the force-extension curve of single stranded RNA. Force- 
extension curves have been measured for single-straned 
RNA and folded RNA molecules [23| but not yet for long 
single-stranded chains. Such measurements have been 
performed for single-stranded DNA strands of about 10 4 
bases 24]. A typical set of results, taken from [24[ is 
shown in Fig. 1131 The force-extension curve rises rapidly 
when the extension reaches the contour length. This is 
preceded by a range of extensions for which the force 
is relatively constant, which may be due to progressive 
loss of stacking interaction. The force-extension curve 
for small extensions is presumably dominate by entropic 
elasticity. We fitted the measured data with a fifth order 
polynomial form, from which we were able to extract the 
quantities v\ (in pN//zm 2 ) and v 2 (in pN//im 4 ), shown 
in Figs. HI and [Til 




1 2 3 4 5 

Extension ( jjuti) 

FIG. 13: Force versus extension data for a segment of single- 
stranded charomid DNA |2fil | with a backbone length of 5.7 
/im in one millimolar phosphate buffer. The data were ex- 
tracted from Figure 1 of |24||. Also shown in this figure is a 
fifth order polynomial fit to the data. 



The V2 control parameter is indeed negative for exten- 
sions less than 3 /im. It should be noted that v\ and v 2 in 
general are expected to be quite sensitive to solvent con- 
ditions. We can now make use of the methods described 
above to compute the Poisson ratio for given ratio I /a, 
and we find that there is indeed a negative Poisson ratio 
over the full range of values of the length parameter o, 
the value oil/ a having been fixed at 0.45 (see Fig. 1161) . 

These results are, of course, merely illustrative, but 
they demonstrate that RNA should be a beautiful "lab- 
oratory system" for the study of gels with anomalous 
elastic properties. As Fig. ED shows, one can tune the 
control parameter v 2 by adjusting the extension, which 
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2 3 4 5 

Extension (ym) 

FIG. 14: The quantity v\ as defined in 12.51 and as derived 
from the fitting curve in Fig. 1131 




3 4 5 

Extension (pun) 



FIG. 15: The quantity as defined in 12. "il and as derived 
from the fitting curve in Fig. 1131 



can, in turn, be achieved via changes in the osmotic 
swelling pressure of the system. The second key param- 
eter, the l/a ratio, can be "programmed" into the RNA 
molecules by alternating random sequences with com- 
plementary homopolymer sequences of prescribed length 
(e.g. strings of C monomers alternating with strings of 
G monomers). The most challenging feature would be 



a ( \xm) 




FIG. 16: The Poisson ratio, /i, derived from the parameters 
v\ and V2 displayed in Figs. I14l and [l31 The ratio l/a is fixed 
at the value 0.45. Note that the parameter a on the horizontal 
axis of this plot refers to the distance between rod centers and 
does not coincide with the extension shown in Figs. 1141 and 

on 



to generate alignment between the rods to produce a ne- 
matic elastomer. Curiously, large ss RNA molecule n the 
form of 1,400 base- long viral genomes are found to be 
quite anisotropic according to low angle X-ray diffraction 
and light scattering studies j2|| . This suggests that large 
ss RNA molecules may be naturally anisotropic. Assem- 
bly of the gel under mild shear flow may enhance this 
natural anisotropy. In summary, RNA gels are expected 
to be rich laboratory systems for the study of fundamen- 
tal elasticity 



A. Implications for viral assembly 

It is interesting to speculate about possible conse- 
quences of a negative Poisson ratio for RNA gels. A curi- 
ous problem of viral assembly is the large discrepancy be- 
tween the density of the genome in solution and the same 
genome encapsidated by the virus's protein shell. Encap- 
sidation proceeds by self-assembly: a solution of capsid 
proteins and viral RNA molecules will spontaneously as- 
semble to form infectious viruses under physiological con- 
ditions. The assembly is driven by generic electrostatic 
affinity between the proteins and the RNA molecules, 
though specific interactions are required to initiate as- 
sembly see, for instance 27] . The scenario for the encap- 
sidation of single-stranded RNA genomes into spherical 
viruses is not fully understood, but Fig. El indicates 
a likely scenario. A partially condensed RNA molecule 




FIG. 17: The first steps in the assembly of a spherical RNA 
virus. In Fig. 117b . the capsomeres are beginning to come 
together. They attach to the portion of the genome in imme- 
diate proximity to them. Figure IT7b shows what happens to 
the RNA complex as the capsomeres assemble and the RNA is 
compressed horizontally if the elastic response of the complex 
is characterized by a positive Poisson ratio. Figure [T7fc shows 
what happens if the RNA complex has a negative Poisson 
ratio. 



forms a condensation surface for capsid proteins of op- 
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posite (i.e. positive) charge, and a curved shell starts 
to form. The challenge for the assembly process is how 
compaction of the genome can be achieved. If the RNA 
material had a positive Poisson ratio, then compaction of 
the genome at the nucleo-protein interface by the electro- 
static affinity would produce swelling of the remainder of 
the molecule that is not in close contact with the grow- 
ing shell. It would seem that completion of assembly is 
not possible. A negative Poisson ratio, on the other hand, 
would lead to collapse of the RNA material into the shell. 
It would be fascinating if model systems could be devel- 
oped with, for instance, RNA gels in contact with a pos- 



itively charged substrate to verify this scenario. We con- 
clude that RNA gels with negative Poisson ratios would 
be well-adapted for easy encapsidation. 
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